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There is convincing numerical evidence that fractional quantum Hall (FQH)-like ground states 
arise in fractionally filled Chern bands (FCB). Here we show that the Hamiltonian theory of Com- 
posite Fermions (CF) can be as useful in describing the FCB as it was in describing the FQHE 
in the continuum. We are able to introduce CFs into the FCB problem even though there is no 
external magnetic field by following a two-stage process. First we construct an algebraically exact 
mapping which expresses the electron density projected to the Chern band, Pfcb, as a sum of Girvin- 
MacDonald-Platzman density operators, pgmp, that obey the Magnetic Translation Algebra. Next, 
following our Hamiltonian treatment of the FQH problem, we rewrite the CMP operators in terms 
of CF variables which reproduce the same algebra. This naturally produces a unique Hartree-Fock 
ground state for the CFs, which can be used as a springboard for computing gaps, response functions, 
temperature-dependent phenomena, and the influence of disorder. We give two concrete examples, 
one of which has no analog in the continuum FQHE with v = | and a xy = |. Our approach can 
be easily extended to fractionally filled, strongly interacting two-dimensional time-reversal-invariant 
topological insulators. 

PACS numbers: 



INTRODUCTION AND STRATEGY 



Models with no net magnetic flux but with a quantized 
Hall conductance a xy have been known since the work of 
Haldane 1 and Volovik 2 . The breaking of time-reversal 
symmetry, necessary for <j xy =^ 0, manifests itself as a 
nontrivial Berry flux for the band, whose non-zero inte- 
gral over the Brillouin zone (BZ) gives the Chern number 
C . The work of Thouless at aP, equates C to the dimen- 
sionless Hall conductance of the filled band. We use a 
convention in which a xy = —C. 

While we focus on single Chern bands, the approach 
to be described here readily applies to strongly interact- 
ing two dimensional time-reversal invariant topological 
insulators (2DTI's^l wn i c h can be thought of as pairs 
of time reversed Chern bands. 

A question that has recently attracted much attention 
is whether these FCB's could also exhibit the FQHE at 
partial filling in the presence of suitable interactions. In 
such cases they are called fractional Chern insulators, or 
FCIs. Optimal conditions call for a hierarchy of scales, 
where the band gap A, the interaction strength V ee , and 
the FCB bandwidth W obey i > y ee » W. There 
have been three fronts of attack. Numerical efforts have 
concentrated on "flattening" the FCEP"^, and realized 
Laughlin-like states by exact diagonalizatiorP^. Most 
recently, other principal FQH fractions such as 2/5 and 
3/7 have been seen as welP^l On the analytical front, 
QfHl has constructed a basis in which known FQHE wave- 
functions can be transcribed into the FCB. Recently Wu, 
Regnault and BernevigS^ have pointed out that if Qi's 
plan is to yield wavefunctions with substantial overlap 
with the exact functions, his proposal must be modified 
to exploit a residual gauge freedom that maximizes the 
overlap. Several studies have likewise been devoted to 
the parton construction for FCl d 16 " 18 * in which the elec- 
tron is fractionalized into quarks, each of which is in an 



Integer Quantum Hall state. 

Our work was stimulated by the third approach due 
to Parameswaran el aP^ who examined the algebra of 
p FCB (q), the density operators projected into the FCB. 
Recall that in the LLL problem the projected density is 
essentially PoMp(q), the Girvin-MacDonald-PlatzmarPSl 
operator, which obeys the algebra of magnetic transla- 
tions: 



[pGMp(q),PGMp(q')l = 2 * sin 

where 

I = 



q x q 



p G Mp(q + q')- 



i 



(1) 

(2) 



is the magnetic length associated with the perpendicu- 
lar external field B . By contrast, the algebra of p F cB(q) 
does not even close, though in the small q, q' limit the 
commutator is proportional to q x q'. The fundamen- 
tal reason for the non-closure of the density algebra is 
the varying Chern density £>(p) in the Brillouin Zone. 
Parameswaran et al offer interesting ways to combat the 
varying B, such as smoothing it out or replacing it by its 
average. 

Our approach, by contrast, tackles the varying Chern 
density from the outset. It is based on two indisputable 
facts: 

• We are looking for the FQHE in the FCB problem. 

• Composite Fermions are very useful in describing 
the FQHE problem in the continuunPH. 

It is then reasonable to ask if CFs can be made play 
an equally fruitful rule in the FCB problem. Our an- 
swer is an emphatic yes. We employ the Hamiltonian 
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approach^!, which provides an operator realization of 
CFs. In the past this has allowed us to compute not 
only gaps but also correlation functions at non-zero q 1 ui 
and T and disorder. At v — | it yielded relaxation rates 
and polarization as a function of T . We describe how 
these ideas can be imported to the FCB problem. 

The Hamiltonian theory of CFs presumes the existence 
of a uniform external magnetic field B . To those who 
say "Where did the magnetic field come from?", we say 
"Where did it go when the v — \ state was described as 
a Fermi sea?"^. We adopt the pragmatic view that one 
must choose whichever mapping takes us closest to the 
desired end product. In the present case, when we want 
to describe FQH-like physics in a Chern band, LL-based 
constructs are a natural platform from which to make the 
leap. Furthermore, as shown by the recent work of Wu, 
Jain, and SuiP^, the Hofstadter problem is adiabatically 
connected to the Chern band problem. 



A. Brief history of the Composite Fermion 

Let us begin by asking what we mean by the CF, since 
the term has many connotations. 

It all began with the realization that in two dimensions 
the statistics of particles could be altered by a singular 
gauge transformation of the wavefunctions that essen- 
tially attached flux tubes to the particles^. For fractions 
of the form v = 2s 1 ) _ 1 , Zhang, Hansson and KivelsorPsl 
converted the electron to a composite boson by attach- 
ing 2s + 1 flux quanta in the path-integral formulation of 
a Chern-Simons theory. They then explained many of the 
FQHE effects in terms of Bosc condensation. JairpBthen 
discovered that when v = „ p . = , one could get excellent 
trial wavefunctions by converting an electron to another 
(Composite) fermion by attaching 2s flux quanta. Lopez 
and Fradkin 27 implemented this flux attachment for Jain 
fractions in a Chern-Simons path integral and computed 
response functions. These flux-attached CFs live in the 
full fermionic Hilbert space, have a bare mass, and carry 
the same charge as the electron, i.e, e* = e. They proved 
especially useful in the gapless state at v = |, analyzed 
in depth by Halperin, Lee and ReacP^. 

The mean-field wavefunction due to flux attachment 
for fractions of the form v — ^r-^-r is 

2p+l 

~^ \z t - Zj\ 2 

where Xp stands for p-filled CF-LLs. Jain's ansatz 

**>*=pJ[(z i -z j ) 2 X p(z,2) (4) 

i<j 

is obtained by dropping \zi — zj\ 2 and projecting the 
z in Xp t° the LLL using V : z — > 2l 2 d/dz. The double 
zero in the analytic Jastrow factor describes the charge 
deficit due to a double vortex that follows the electron, 



leading to a CF that has e* = e(l - = ^pj- This 

is the CF obtained by vortex atta chmen t. 

In the path integral approache d 26 * 27 * the change from 
flux attachment to vortex attachment is achieved by con- 
sidering fluctuations about the mean-field, while in our 
earlier Hamiltonian approach^] plasmon correlations a la 
Bohm-Pines were responsible. 

B. Brief review of the Hamiltonian theory 

We worlP^l with CFs that live in the LLL from the 
beginning, as did ReacP^, and Pasquier and Haldane^. 
Our CFs carry both the phase and charg e deficit of a 
double zero in the FQHE wavef unct ior J 2 1 * 3 M . Their entire 
Hamiltonian is given by the electron-electron interaction 
projected to the LLL. 

We introduce CFs as follows. The FQH problem pro- 
jected to the LLL is defined by the Hamiltonian 

H = W XI ^( ( l) l, ee(q)PLL L (-q) (5) 

q 

where jO LLL (q) is the electron density projected to the 
LLL. In first quantization the full electron density is 

p(q) (6) 

3 

The electron's position r e may be decomposed as 

r e = R e + r] e (7) 

where the electronic guiding center coordinate R e and 
cyclotron coordinate rj e obey 

[Rex,Rey] = —H 2 (8) 
[VescVey] = (9) 

[r, e ,R e ] = 0. (10) 
Upon projecting to the LLL 

Plll (q) = E^-^e^ = e -^ 2/ VGMp(q). 

3 

(U) 

where each term e !q Rcj in the sum obeys the GMP 
algebra by itself thanks to Eq. [8] We shall use the same 
symbol for the densities when we switch to second quan- 
tization. 

So the Hamiltonian to solve is 

h = 5 PGMp(qK e (q)/>GMp(-q) (12) 

q 

where v ee (q) — v ee (q)e~ q2[2 1 2 . The problem is difficult 
because, with r] e projected to the LLL, the electron is 
described by just one canonical pair R e and not two. The 
LLL projected electron has half the degrees of freedom of 
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a regular two-dimensional fermion. However, the biggest 
problem is that at fractional filling there is no clear mean- 
field state. 

We attack these problems as follows. First we intro- 
duce an auxiliary pair of conjugate "vortex" guiding cen- 
ter coordinates R u . They are defined by their commuta- 
tion relations: 



[-^;a; 3 Rvy\ — 



J 2 



2p 



2p+l 



(13) 
(14) 



Evidently the vortex describes a particle whose charge 
— 2 p+i m electronic units is exactly that of the two vor- 
tices in the Jastrow factor. It too has just half the degrees 
of freedom of a regular two-dimensional particle. 

We want these auxiliary coordinates to commute with 
everything electronic i.e., 



[R e , R„] — 0. 



(15) 



The cornerstone of our approach is that we can accom- 
modate both R e and R„ and their algebra very neatly into 
the Hilbert space of a regular two-dimensional fermion, 
which is going to be our composite fermion. This 
fermion feels the reduced field B* seen by a charge e* 
object. In terms of its guiding center and cyclotron coor- 
dinates (which carry no subscripts like e or v) that obey 



[Vx,Vv\ 



il* 



i 2 



1-c 2 



[R x ,R y ] — —il* 
[T7,R] = 



(16) 

(17) 
(18) 



the algebra of the two conjugate pairs R e and R„ can be 
realized as follows: 



How is this formalism, predicated on making the CF- 
substitution in p GMP , to be applied to the Chern band 
problem where the density of interest is /3 F cb? The key 
is to establish the following algebraically exact mapping: 



p F c B (q) 



G 



c(G,q)p GMP (q+ G) 



(22) 



where the coefficients c(G, q) can be computed from the 
data on the original Chern band, essentially by Fourier 
transformation. The CF-substitution can be then made 
in each /9 GMP (q + G). While an explicit demonstration 
follows later, here is the gist of the argument. On an 
N x N toroidal lattice the number of possible values for 
p and q are TV 2 each. We will show that the /9 GM p(q+G) 
are linearly independent only for iV 2 values of G for each 
q restricted to the Brillouin Zone (BZ). These N 4 linearly 
independent operators p GM p(q + G) form a complete ba- 
sis for one-body operators, just like the canonical basis 
dt( Pi )d( Pi ) [i,j = l...N 2 ]. 

The outline of the paper is as follows. In Section II we 
will show that our approach applies to a (Type I) Chern 
band with a variable Chern density B(p) and coulomb 
interaction between electrons. This Chern band is ob- 
tained by starting with two electronic LLs and applying 
a periodic potential V PP which mixes the LLs and causes 
£>(p) to vary. The lower band, the Modified Lowest Lan- 
dau Level, or MLLL is our Chern band with C = — 1. 
We explicitly derive Eq. 22 (with the MLLL being the 



FCB), compute the coefficients c(G,q) and calculate the 
band structure of CFs in the HF approximation. 

In Section III we show how, given a specific lattice 
Chern band and an interaction between electrons, one 
can introduce CFs without any reference to LLs. We 
choose as our Type II example the Lattice Dirac Model 
(LDM): 



R e 
R„ 



R + 77 c 
R + rj/c. 



(19) 
(20) 



This in turn permits the crucial CF substitution 

PcMp(q) = X] exp t iq ' ( R -? + cr? ^] ( 21 ) 



in Eqn. 
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for the projected Hamiltonian H , which 
now acts on a regular fermionic Hilbert space with two 
conjugate pairs per particle. Since the CFs see exactly 
the right field to fill an integer number of CF-LLs, a 
natural, gapped Hartee-Fock state emerges. The price 
we pay for obtaining a good mean-field starting point is 
that our Hilbert space has unphysical degrees of freedom 
R„ . In order to work in the physical sector the vortex co- 
ordinates need to be constrained. Specifically, the vortex 
densities, p v (q) — e lq R " , emerge as a gauge algebra. The 
way to handle this gauge degree of freedom is described 
in our review^. 

The numbers computed in this scheme at nonzero 
ui,q,T and disorder agree with data at the 10-15 % 
levePl. 



H(p) = a 1 smp x +a2 sm.py-\-a^(M— cosp x — cosp y ) (23) 

with M = 1, which lies in the regime with C = — 1. 

Models in which LLs appear or do not appear explic- 
itly are labeled as Class I and Class II respectively. 

In Section IV we ask what happens if we apply our 
approach to a band with C = 0, while in Section V we 
turn our attention to FQH-like states that owe their very 
existence to the lattice potential. Flux attachment on 
a lattice was first investigated in the context of anyonic 
statep2122l and analyzed in the FQH context by Kol and 
ReacP^. In such states, due to the explicit breaking of 
Galilean symmetry, the dimensionless Hall conductance 
need not be equal to the filling factor. There is suggestive 
numerical evidence of such states in a p roble m of hard- 
core bosons in an external magnetic 

fiel( j35|36| By yirtue 

of our mapping of the FCB problems to LLL problems, 
such states should exist in FCBs as well. We define and 
solve an illustrative Class I model at v = \ and demon- 
strate the existence of such a state where v = ~ and 

5 

_ 2 

u X y — 5 ■ 

Section VI presents conclusions and open questions. 
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II. CLASS I MODELS 

The goal of this section is to convince the reader that 
a nonconstant B(p) is no impediment to the CF substi- 
tution, and to flesh out the key expansion Eq. [22] We 
begin with the construction of a nontrivial Chern band 
with a non-constant £>(p). 




MLLL 



The simultaneous eigenfunctions we seek are^: 

(x e , y e \p, n) = ^ p , n (x e ,y e ) (27) 
V jy-toyWle^Mxt-aj-pyl 2 ) (28) 



j=-oo 



where <p n (x e — aj — p v l ) is the wavefunction for an os- 
cillator in level n centered at x e = aj + l 2 p y . 
Hereafter we will set a = 1 which means 



I 2 



1 

2V 



(29) 



The states are normalized to unity, and r e integrals go 
over the unit cell. 

The Bloch functions are 



Kp,n)) = e-^\p,n) 
and the Berry connection 



(30) 



FIG. 1: Left: Two unperturbed Landau Levels with C = 
— 1 and zero width separated by energy lu. The dotted grid 
represents a fictitious square lattice with one flux quantum 
per unit cell. Right: The two bands with finite width after a 
periodic potential V PP (solid grid) is imposed. The lower of 
the two bands, the Modified Landau Level, the MLLL, is our 
Chern band with C = — 1. 



Consider a problem with two LLs labeled and 1 sep- 
arated by a gap oj that is at our disposal, as shown in 
FigjT] By choosing the Hamiltonian to be a suitable func- 
tion of r]\ri e , not simply linear, we can arrange for the 
other LLs to be separated by a parametrically larger gap 
than uj and hence ignorable in what follows. 

Each level has C = — 1. It is instructive to demon- 
strate this explicitly. (We recommend the review by 
Xiao, Chang and Niu for some basic ideas of magnetic 
Bloch bands^l.) First we mentally superpose on this con- 
tinuum immersed in a perpendicular field Bq, a square 
lattice of side a. No real periodic potential is applied yet. 
Working in the Landau gauge 

Ay(x,y)=xB Q A x (x,y)=0 (24) 

we seek energy eigenfunctions which are also simultane- 
ous eigenfunctions of T x and T y , the magnetic translation 
operators in the x and y directions: 



(25) 



These commute with H, but not with each other unless 
each unit cell has an integer number of flux quanta. We 
choose the simplest case of one flux quantum penetrating 
each unit cell, i.e., 



a 2 = 2ttZ 2 . 



(26) 



A(p ,n) = i(u(p,n)\V p \u(p,n)) (31) 
can be computed to have components 



A y = A x = Pyl 2 = —p.y so that (32) 
B(p) = V p x A = — — which means (33) 

ZTT 



2tt 



B = -1. 



(34) 



BZ 



However this B is constant in p in both LLs. To make 
it vary, we add a periodic potential 



V(r e ) 



G 



V(G)e 



G 



(35) 



which mixes the LLs and induces structure in B(p). 
In our illustrative example we keep only the harmonics 
±27r in the two directions with coefficient V\o, though 
the following analysis applies to the general case. Using 



(pn 2 \e lG - r \p ni ) = Pn2ni e^ G y/^e tGx ^ (36) 



where (for the general value of I), 



y n >- 

\ |m-n 2 | 

t|J when n x > n 2 

/ , \ lm-«2| 

1 73/ when 712 ~ Ul 

z = q x + iq y 



q 2 l 2 
2 



(37) 
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we find 



H T (P)= 9(p) - Vai sin p y + Va 2 sin p x 



/ttV 



(cos p x 



iPy) 



V 



V 10 e~ 



(38) 
(39) 



The ground state of H 1 (p) will be referred to as the 
Modified LLL, or MLLL. ft is our Chern band. The 
function g(p) affects the energy dispersion of the MLLL, 
but not the Chern density B(p). 

Though H 1 (p) has the form of the Lattice Dirac Model 
(Eq. 23 1 , it is in the topologically trivial region. This is 



due to our requirement u) > 1~J~kV which ensures that 
the two bands do not touch at p x — p y = 0, which in 
turn ensures that the Chern number remains C = — 1. 
Due to the topological triviality of i? 7 (p) the pseudo- 
spin n(p) = (u(p)|<x|u(p)) never enters the southern 
hemisphere. Nonetheless the overall C = — 1 because 
nontrivial topology is contained in the p-dependent ba- 
sis functions. Whereas in the traditional LDM, the 
tight binding wavefunctions arc p-independent spinors, 
[1,0] T and [0, 1] T , here they are the states |p, n = 0, 1) 
with topologically nontrivial p dependence. The total B 
in this problem has a constant piece — ^ coming from 
the basis functions and responsible for C = — 1, and two 
more p-dependent terms with zero integrals: one due to 
the p-dependence of the ground state spinor, and a cross 
term that arises because (n|VJn') 7^ for n 7^ n' . The 
total B(p) is shown in Figure [2] along with the Chern 
density for the LDM at M = 1. Notice the strong simi- 
larity even in this minimal model with just two LLs and 
one harmonic in V. 

Now that we have a nontrivial Chern band with a 
nonconstant B(p) let us proceed to the CF-substitution, 
which will in turn lead us to the gapped state in the HF 
approximation at v = | when interactions are turned on. 

First we need to find Pmlll, the projection of the elec- 
tron density operator to the Chern band, the MLLL. 
When V10 = 0, clearly p ML L L (q) = Pux(q)- To find out 
what it is when V10 is turned on we proceed as follows: 

• Find the 2x2 matrix that describes the electron 
density p e (q) = e iq ' rc in the space of the LLs, n = 
0,1. 

• Find the eigenstates of H I (p). 

• Project yO e (q) = e lq l " e to the ground state at each 
p, the MLLL. 



The matrix elements of /9 e (q) 



between the 



magnetic Bloch states defined in Eq. [28] vanish unless 
the initial momentum p and final momentum p', both 
restricted to the BZ, obey 




FIG. 2: Top: B(p) in the MLLL with just two LLs and one 
harmonic Vol 7^ 0. Bottom: Berry flux density in the Lattice 
Dirac Model at M = 1 . 



2wN x e x + 2nN y e y = p + q [p + q] (41) 
The non-zero matrix elements are found to be 



([p + q]n 2 | e ^ r =|pn 1 ) = 

\ \lxly + QxPy - (Px + QxfaNy 



Pn^rix CXp 

EE p Il2ni (q)e l *^P) 



(42) 



where p n2ni (q) has been defined in Eq. 37 The asym- 
metry between p x and p y in Eq. 42 reflects our choice of 
the Landau gauge in defining the basis states in Eq. |28| 

In view of its importance to what follows we display 
e »*(q, p) prominently below: 



p' = [p + q] = (p + q) mod G 



(40) 



Thus we must subtract from p + q that G which re- 
stricts p' to the BZ. 



e **(q. p) = ex p 



l 

2^K 



:Q_ x q y + q x p y - {p x + q x )2nN y 



(43) 
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The corresponding second-quantized operator is 



P ni, 7i2— 0,1 

(44) 

where a n and a' n are the operators associated with the 
basis states |pn). 

Let U be the matrix that diagonalizes -ff 7 (p) in Eq. 
[38] and relates a n to the d n associated with the energy 
eigenstates as follows 



«o(p) 
ai(p) 



C/io 



C/oi 



(45) 



Since H 1 (p) is topologically trivial, U, like the eigen- 
spinors, is globally defined in the BZ and periodic in p. 
Switching to the new basis and projecting to the ground 
state we obtain 



p MLLL (q) = E4(p')4(p)e 4 * (q ' p) /(q,p) 
p 

/(q,p) = c4,(p'K'„(q)^o(p) (46) 

Hereafter the subscript on do, indicating that it corre- 
sponds to the ground state or MLLL, will be dropped. 

Thus we have a Chern band, a non-constant B and 
a closed expression for the projected density. The final 
step before we carry out the CF substitution is to write 
this density in terms of p G mp- Before doing this explic- 
itly, we provide an intuitive argument that this can be 
done. When V w = 0, we know p M lll = e~ q 1 'Vgmp- 
As we turn on Vio, the perturbing terms are of the form 

„iG -r e 



e iG -T7 ee iGo-R e where Q Q _ 2n (e x 



lyUy) 



with only one of n x or n y = ±1. Given the GMP algebra, 
the repeated action of this perturbation can only be to 
turn pcMp(q) into a sum over p G Mp(q + G), where G is 
now any reciprocal lattice vector. So we do expect that 
in the end, even for an arbitrary periodic potential 



Xq) = E c ( G ' q )'° GMp ( q + G ) 



(47) 



G 



Since the bands never touch, perturbation theory 
will always converge. However the final result is non- 
perturbative and follows simply from the dependence of 
ff 7 (p) on e iG °- R -. 



We will now show Eq. 47 explicitly and compute 
c(G,q). 

First let us construct an auxiliary operator which 
obeys the GMP algebra. 



PoMp(q) = E d V) d (P) el * (q ' P) 



(48) 



Given any BZ in which the operators d, d^ appearing 



in Eq. 48 are canonical, it is easily verified that this 
operator satisfies the magnetic translation algebra Eq. 
[I] We likewise construct operators with momenta 



q + G = q + 2im x e x + 2irn y e y : (49) 
defined by 

p G Mp(q + G) 

= £V(p / M(p)e i * (q+G ' p) 
p 

_V^t (p / )(i(p)e i *^ q ' P- 1 e i(ly n *-1* n y+ 27rn * n y) e -iPxn y +ip y n x 

(50) 

Note that whether we transfer momentum q or q + G 
to p the resultant p' is the same. We emphasize that 
these operators can be constructed for any problem in a 
square lattice BZ, with no reference to any LLs. This fact 
will be crucial in the next section. 

We now give the details of the counting argument that 
assures us that /9 MLLL (q) may be expanded in terms of 
PcMp(q + G). Consider a system of size L x L wrapped 
into a torus. The question to ask is: For a given q in the 
BZ, for how many different values of G are the /0 G Mp(q + 
G) linearly independent? Since a = 1, the number of 

sites is TV 2 = \ — L 2 , which also equals the number 
of points in the BZ, the number of distinct values for 
p and the number of distinct values of q in the lattice 
model. The smallest value for any component of q or p 
is q m in — Pmin — • To verify that the largest distinct 
value for any component of G is G max = 2-kN consider 
the second and third exponentials in Eq. [50] which alone 
depend on G. Focus on a factor like e~2 9=™B when q x 



ami n„ 



e 2 



N 



(51) 



The same goes for all the terms in the second exponential, 
while the third exponential always equals unity, which 
means 



Pgmp (q + G max ) oc pcMp(q) 



(52) 



Thus we get linearly independent densities only for 
components up to G max = 2nN. There are only iV 2 
independent values of G, just as for p or q. But this 
means there are iV 4 linearly independent operators of 
the form /? G Mp(q + G), exactly the right number to 
form a basis, like the canonical basis d^ p2 d pi . So what 
we find is that not only p FC B, but any bilinear oper- 
ator O (such as the lattice current operator) of the 
form ^ p eft (p')d(p)0(q, p) can be expanded in terms of 

p G Mp(q + G). 

Having hammered home our central point, we now 
turn to the determination of the coefficients of the ex- 
pansion. To this end we combine Eqs. [46} [47| and [50| and 
write 



7 



The CF substitution 



PMLL L (q) 

= 5^ rft (P , )d(p)e i * (q,p) /(q,p) 
p 

= E E c(n x ,n v ,q)dt(p / )d(p)e < * (q ' p) 



This equation can of course be satisfied since 



(53) 



(54) 



/(q,P) = E c ( n x,n y ,q) 



is, at each q, just the Fourier expansion of the function 
/ periodic in p in terms of oscillating exponentials of the 
right period. 

The commutators of the projected electron density 
(q) can be worked out, if desired. They will be 
neither pretty nor universal, unlike the magnetic trans- 
lation algebra^, depending instead on the details of the 
lattice via /(q, p). 

Having expressed /0 M lll(ci) in terms of p G Mp(q + G) 
we need to do the same for a term in H that is absent 
in the usual LLL: the non-constant kinetic energy — e(p) 
of the MLLL. Because e(p) is periodic mod G, this is a 
special case (q = 0) of the Fourier transform we carried 
out for jOmlll (q) ■ We write 



5>t(p)d(p) £ (p) 
p 

= ]TMGKmp(G) (56) 

G 

= d^p)d(p)h(n x ,n y )e- tp * n y +ip » n * +27rn * n i57) 



which amounts to Fourier expanding the energy dis- 
persion — e(p). We now have the full electronic hamilto- 
nian for the FCB in terms of p GM pS- 

Note that the phase factor 



o**(q. p) 



exp 



1 

2tt V 2 



q x q y + QxPy - ip x + q x )2nN i 



(58) 

jumps in p space for a fixed q. For example if q = 3e y 
and the BZ is in the interval [0 — 2tt] in both directions, 
then for any p with p y > 2tt — 3, adding q will take it 
to the next BZ and N y will have to jump from to 1 
to bring p' within the BZ. Luckily, this discontinuous $ 
and its jump are shared by both /0 MiM (q) and /9 OM p(q + 
G) as we find in Eqs. 53 and 54 This ensures rapid 



convergence of the Fourier expansion of the jump-free 
part /(q, p) in Eq. 55 



Now we must switch to CFs. We only sketch the broad 
ideas. We consider the case of v = | when the CF has a 



charge e* 



|e and I* 2 



3l 2 



: . The spatial unit cell 



is 3 units long in the x-direction so as to enclose unit flux 
as seen by the CF so that BZ CF goes from — ? < p x < | 
and is unchanged in the y-direction. However when we 
construct the projected electron density operators we will 
need to consider q that runs over the BZ of the electron 
not the CF. 

Consider p GM p(q) which was e iq Rc in first quantiza- 
tion and 



/OcMp(q) 
= <*VMp) 



p£BZ c 



x exp 



1 

2tt V2 



q x q y + q x p y - (Px + q x )2irN y 



(59) 



in second quantization. To go to the CF representation 
means to write 



iq-R e _ iq-(R+T)c) 



(60) 



in first quantization and the following representation 
in terms of CF operators C and in second quantization 



PcMp(q) 

= E ^,(p / )C n (pK/ n (q->oq ) /->r) 



PGBZcf 

3i /l 



x exp 



2tt V 2 



q x q y + q x p y - (Px + q x )2-KN y 



(61) 



where the 3 is due to I* 2 = 3/ 2 = and the argument 
of Pn'n is qc because the c in r\c may be lumped with q 
(see Eq. 21 ). Note that all CF-LLs (n — 0, 1, . . . ) appear 



in the density, a result of the enlarged Hilbert space in 
which we are representing the problem. 

With the Hamiltonian expressed in terms of CF op- 
erators, we move to the HF calculation. We set uj = 10 
(the gap between the two electronic LLs ) and choose 
the periodic potential to be V PP — Vioe~^' K — 1. We 
keep 3 CF-LLs which get mixed by the periodic poten- 
tial and interaction. We end up with three bands which 
are fairly well separated. Unlike in the continuum where 
the CF-LLs were uniformly filled, the occupation num- 
ber here varies with p (due to the periodic potential) 
and had to be found self-consistently. Figure [3] shows the 
results of our calculation for the coulomb interaction of 



strength 



2ttV c , 



We see a clear gap separating the low- 



est band which is fully occupied from the others, even at 
very small values of V ee . 

This is puzzling since one expects the CF-picture to 
break down for V ee <C V rP . Upon further investigation 
we found two signals that point to the breakdown of the 
CF picture, one internal and one external. 



Hartree Fock Bands, nu=l/3 




FIG. 3: The results of our Hartree Fock calculation at v = | 
and lu — 10, for the Coulomb interaction 27r ^" . The three 
bands resulting from three CF-LLs which get mixed and mod- 
ulated by the periodic potential Vppp = 1. 

The internal one involves the occupation numbers of 
the CF-LLs in the ground state at each p. Figure [4] shows 
that at V ee = 10, uqf = 0, 1 are robustly occupied while 
Tic f — 2 has negligible occupancy. Thus our truncation 
with three CF-LLs is safe, since the n — 2 level is not 
called into play. By contrast at V ee = 0.5 we see in Figure 
[5] that occupation of Ucf — 2 can be substantial. The 
levitation of the fermions to the upper CF-LLs is a clear 
indication that CF theory is failing as a good low-energy 
theory. 

Before moving to the external signature, let us observe 
a three-fold symmetry in the occupations as a function 
of p y , which is also reflected in the energy bands. This 
is a consequence of the x-translation symmetry of the 
nonintcracting HamiltoniarP3. Recall that T x commutes 
with the Hamiltonian, but not with T y . It can be easily 
shown that 

T x \ Px , Py ,n) = e^\p x ,p y ~Y' n ) ( 62 ) 

This symmetry is also possessed by the HF Hamilto- 
nian 

T^.H HF (p x ,p y )T x = H HF {p x ,Py - y) (63) 

and thus by the energy bands and the occupations. 

Now we turn to the external test that signals the 
breakdown of the CF state and is sharper than the one 
based on occupation numbers. It involves the compari- 
son of the variational energy per particle in the FCI state 
and the electronic Fermi liquid state. It is clear that the 
energy of the elecronic Fermi liquid state (which lies in 
the physical Hilbert space) provides a variational upper 
bound for the exact ground state energy. It is not so 
clear that the energy of the CF state, which is defined 
in a bigger space containing the physical coordinate R e 
and the unphysical vortex coordinates R„ , is variational. 
One might worry that these extra unphysical degrees of 
freedom may allow a further lowering of energy not per- 
mitted in the physical sector. This is, however, not true. 
Since the hamiltonian in the enlarged space depends only 



Occupations n=0,1.2. Vee= 10.0, Vpp=1 




FIG. 4: The occupation numbers in CF-LLs 0,1, and 2 at 
Vee = 10,Vp P = 1, uj = 10. Notice that the n = 0, 1 levels 
(bottom two) saturate the occupancy while the n = 2 (top) 
level is practically empty, validating the truncation at three 
CF-LLs and establishing CF theory as a good low-energy de- 
scription. The plots repeat three times in the p y direction 
because T x and T% commute with H, but do not commute 
with T y and add ^ and 4^ to p y , respectively. 



on the electronic guiding center coordinate R e via p G mp 
and is independent of the vortex coordinate R„ , the ex- 
act eigenfunctions must be tensor products of the exact 
eigenfunctions in the physical sector and arbitrary wave- 
functions in the unphysical sector. However, the energy 
is independent of the choice made in the unphysical sec- 
tor. Thus the exact ground state energy in the enlarged 
space is equal to that in the physical sector and conse- 
quently any function in the enlarged space can furnish 
a variational upper bound to the exact ground state en- 
ergy. We compute the energy of the HF ground state 
with one filled CF-LL. This, of course, is not an exact 
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HF energies: FCI versus Fermi liquid 



Occupations n=0,1,2, Vee - 0.5, VPP=1 




FIG. 5: The occupation numbers in CF-LLs 0,1, and 2 at 
Vee = 0.5, Vpp = 1, u! = 10. Notice that the particles are 
levitating towards the n — 2 level (top) and moving away 
from n = 0, which calls the truncation at three CF-LLs into 
question. The plots repeat three times in the p y direction be- 
cause T x and which do not commute with T y add ^ and 
4p to p v , respectively, commute with H. 



eigenfunction. In fact, it is not even in the form of a 
tensor product between the physical and unphysical sec- 
tors, being instead a linear combination of tensor product 
states. However, by the above argument, it nevertheless 
provides a variational upper bound on the exact ground 
state energy. Fig. [6] shows the the HF energy per particle 
of the Fermi liquid (smaller dots) versus to | FCI state 
(larger dots). We see that the Fermi liquid yields to the 
FCI state at V ee ~ 2.5. We caution the reader that this 



Energy 
particle 
-1.0 

-1.5 
-2.0 
-2.5 
-3.0 
-3.5 



V 



FIG. 6: Comparison of the HF energy per particle of the 
electronic Fermi liquid (smaller dots) versus the | FCI state 
(larger dots). We see that the Fermi liquid, which wins at 
small interaction strength Vee yields to the FCI state around 

Vee = 2.5. 



does not mean that the FCI state unequivocally wins. It 
is possible that there are correlated electronic states with 
an even lower energy that our CF-state. 

To summarize, we produced a non-trivial Chern band 
by transferring the topology to the basis functions in p 
space. These functions arose from two electronic LLs 
mixed by a periodic potential. The projected charge den- 
sity at momentum q was then written as a computable 
sum over G of GMP densities at q + G. The CMP den- 
sities were treated in the Hamiltonian method by the 
replacement R e = R + rjc. Finally a HF calculation 
was carried out using a gapped CF ground state. We 
see that although a FCI state always exists, even for very 
weak interactions V ee , the occupations of the higher CF- 
LLs becomes smaller with increasing V ee , providing an 
internal signal of its stability. A sharper limit for the 
goodness of the FCI state is provided by the comparison 
to the variational energy of the electronic Fermi liquid 
state, which wins for V ee < 2.5, but gets bested by the 
FCI state for larger V ee . 



B. Variants and limitations of Type I models 

A natural extension of the above example is to impose 
more complicated periodic potentials toget more compli- 
cated B(p)s. Earlier in our exploration^ we pursued this 
line of thought so as to reproduce the B(p) of some spe- 
cific lattice model, say the LDM. (We see in Fig. [2] that 
the B(p) of the MLLL with just one harmonic is already 
not too different from that of the LDM.) Our motivation 
was as follows. Let us take the view that the FCB problem 
is defined by (i) B(p) and (ii) the interaction written in 
terms of the projected density p FCB (q). The logic behind 
(i) is that the projected electron coordinate 



R 



FCB\ 



'PP' 



._d_ 
dpfj 



■A*(p)K(p-pO 



(64) 
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has commutators denned by B: 



(65) 



It follows that if we can construct a surrogate band 
with the same B{p) as in a given FCB, then any function 
of the projected electron coordinates will be algebraically 
the same in the two problems. 

Let us review where this line of thought leads in the 
continuum FQH problem. There we represent the elec- 
tron's projected coordinate R e in terms of CF variables 
R and rj. Getting the algebra of R e correctly also means 
getting the algebra of /0 LLL , the projected density right, 
up to a known prefactor. Thus, if (• ■ ■ ) denotes averages 
in the target band, which here is the LLL, 



r e = R e + r) e which implies 

(r e ) = R e while 

iq-r e \ _ ^ e iq-(R c +T7 c )^ 

= (e^^e 1 ^' 



(66) 
(67) 
(68) 
(69) 
(70) 



In other words, the projection of the exponential of r e 
is, up to a known prefactor e~ q 1 / 4 , the same as the ex- 
ponential of the projection because rj e and R e commute. 
This is why if the projected coordinate is faithfully rep- 
resented, so is the projected density. We go over these 
well known facts to highlight the unusual simplicity of 
projecting to the LLL. 

Unfortunately, in the generic FCB problem this is no 
longer true. Let us define two different projected densi- 
ties in the FCB. One is the usual one: 



p FCB (q) = (FCBle^lFCB) 



(71) 



This is the projected density which enters the interacting 
Hamiltonian in the FCB. 

The other is the analogue of the guiding center density: 



p - FCB(q) = e^ CB 
In the FCB problem 

Pfcb (q) ^ C(q)p FCB VHJ 



(72) 



(73) 



because the "guiding center" coordinates R FCB do 
not commute with the analogue of the cyclotron coor- 
dinates. Thus, B(p) is not enough to specify the inter- 
acting Hamiltonian in the FCB completely. One needs 
the expression for /0 F c£i(q) as well. 

However, to lowest order in q and q' we have, in first 
quantization and in p-space, 

[p FCB (q), PFc B (q')] = -« [q X q'] S(p)+higher order terms 

rn (74) 
as pointed out by Parameswaran at ap3. 



So if the Chern flux density B(p) of the surrogate band 
matches that of the lattice FCB, p FCB (q) and p FCB (q) will 
bear a close resemblance to each other but not be equal 
in all respects. Since these are both models anyway, one 
may argue that it is sufficient to get a surrogate that ap- 
proximates the original lattice model and is yet amenable 
to analytic treatment. However we will not pursue this 
approach further since there is a more direct way to ob- 
tain /0 FCB (q) in terms of /0 GMP (q) for an arbitrary lattice 
model, which we now describe. 



III. TYPE II MODELS 

In the previous section we showed how to carry out the 
CF substitution in LL-based models with nonconstant 
B(p). In this section we present a more general approach 
for incorporating CFs in the solution of any given lattice 
model with a Chern band, without any reference to LLs. 
We illustrate this approach by considering the LDM with 
some interaction V ee . 

Here are the concrete set of steps we follow: 

• Construct jO FCB (q), the projected density operator 
in terms of the eigenfunctions of noninteracting 
LDM hamiltonian H(p) and the operators d and 
d> associated with the Chern band at each p. 

• Construct operators obeying the algebra of 
PcMp(q + G) using d and d) as per Eq. 



50 As 



noted immediately after that equation, this can al- 
ways be done. To get the best possible results this 
must be done in A y = gauge. 

• Expand as before, using the complete set of iV 4 
operators p GM p(q + G): 



PFc B (q) = ^c(G,q)p GMP (q+ G). 



(75) 



G 



This will be just a Fourier expansion in 

• Carry out the CF substitution in p GMP and go on 
to the HF approximation. 

We illustrate the above steps with the LDM at M = 1: 
H(p) — ox sinp x + o~2 sin p y + a ^(l — cosp x — cosp y ) (76) 
with ground-state energy 



- e(p) = -^1 + 2(1 - cosp K )(l - cosp y ). (77) 
The generic formula for the projected charge density 



5 (q) = X) rf+ ([P + Qj)rf(p)<[p + q] |p> 



(78) 



where d) , d create and destroy the ground state and 
|p) is the corresponding Bloch spinor. However |p) is not 
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globally definecP over the BZ because C = — 1. Here are 
two choices that work in two different patches: 



where 



|P)? = 
|P>^ = 
cos#(p) 



Sill 



Hp) 



. cos m e im 



1 — cos p.,. — COSPy 
sinj3 x + i sinpj, 



sin 2 + sin 2 p y 



(79) 
(80) 

(81) 
(82) 



The superscript on the kets reminds us that these 
are going to be transformed to a more appropriate gauge 
later. 

The angle </>(p) is ill defined at 9(p) — 0,7r. At 
9(p) = 0, choice |p) 2 is good because the component 
sin ?^-e~^^ vanishes. In a patch where 9(p) = n, the 
choice |p)i is good. Figure [7] shows that there are four 
trouble spots: (0, 0) where the spinor is at the south pole 
9 = it, and points (0, 7r), (it, tt) and (tt, 0) where it is at the 
north pole 9 = 0. We pick the BZ in the range [-§ , ^] 
so that the trouble spots are not at the edges. Patch 1 
is indicated by the solid right triangle and patch 2 is the 
rest. 

The corresponding Berry connections are 



In this final gauge, not only is A y = 0, the two expres- 
sions above, |p)i and |p)2, merge seamlessly along the 
sloped line p y = j — separating the patches and on 
the vertical boundaries of the BZ, which may be glued 
to form a cylinder. However between the lines p y = — f 
and p y = there is a phase mismatch so that we cannot 
roll the cylindrical BZ to a torus without a discontinuity 
at the seam. This had to be so, for otherwise we would 
have a fully periodic Bloch function and C would vanish 
by Stokes' TheorerrP. 

A salient feature of this gauge is that /9 FC b(cO m Eq. 
[78] also has jumps in the sum over p exactly where the 
PcMp(q + G)'s do, rendering this the optimal gauge for 
rapid convergence of the Fourier expansion. 

Figure [8] gives a taste of how the Fourier expansion 
works at a generic value of q = 3e x +3e y . It compares, at 
fixed p y = — ? + 5, the imaginary part of the coefficient of 

^([P + q])^(p) hi PFCB(q) to the approximation in which 
the Fourier sum over G is truncated after 50 harmonics 
(from —25 to 25) in the p y direction and 20 in the p x 
direction. Note that even though q has components in 
both directions, the jump occurs only in the p y direction 
due to the periodicity in the p x direction. 

Figure [9] shows the full landscape of the the real part 
of the coefficient of d'([p + q])d(p) in p F cB(q) versus the 
approximation in which 50 harmonics (from —25 to 25) 
are kept in the p y direction and 20 in the p x direction. 

Having expressed everything in terms of p GUP (q), the 
the CF substitution and HF analysis can be carried out 
just as before and we do not discuss it further. 

A central message of this work is that since in a prob- 



-A?(P) = 4 (i 



cos0(p))V<j 



A{V) = -(l-cos0(p))V^ 



(83) 
(84) 



Once again the superscript on the »4 's above signals 
that this is not yet the final gauge. 

If we imagine the edges of the BZ parallel to e y glued 
together to form a cylinder and then the top and bottom 
sewn together to form the torus, this dark line will mark 
the boundary between the two regions. The difference 
between the two connections is — V0, and the integral 
around the boundary of — ^Vcf> will yield C = — 1. 

We now reach the final gauge A y = as follows. Let 
us define 



^l(Px,P y ) 

A 2 (Pz,P a ) 



■Aly(Px,P y ')dp y 



A i(Px, 4 - y) + 



x(px) = Hp* 



TV 

4 



Px, 



(85) 

■A y i.Px,Py')dpy 

(86) 
(87) 



and the following final kets in the two patches 



\V)x 

|P>2 



e lAl |p)i 
e lA2+lx |p>2 



(89) 



(-jt/2,3Jt/2) 



(-Jt/2 



(-jt/2,-jt/2) 



(3jt/2,3jt/2) 




(3jt/2,-jt/2) 



FIG. 7: In patch 1, we have the south pole at (0, 0) and spinor 
|p)i is well defined while in patch 2, spinor |p} 2 is well defined 
at the north pole reached at the points (0, tt), (tt, 0), (71", 7r). 
The patches meet on the right triangle. In the final A y = 
gauge the wavefunction is seamless across the hypotenuse 
p y — j — ^f- and periodic in p x . However the top and bottom 
edges differ by a phase and the Chern number is resident in 
that difference. 
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FIG. 8: Here we compare, at fixed p x — — ^+5, the imaginary 
part of the coefficient of d) ([p + q])d(p) in pFCB(q) to the 
approximation in which the Fourier sum over G is truncated 
after 50 harmonics in the p y direction (from —25 to 25) and 
20 in the p x direction. 



lem with C ^ 0, one cannot work with periodic Bloch 
functions, the expression for p FCB (q) in the A y — gauge 
will necessarily have a jump beyond some p y depending 
on q y when we retract from a point p + q which lies out- 
side the BZ to a point p' = [p + q] within. The GMP 
density PoMp(q), has exactly such a jump and is the right 
basis to use. When the Chern band supports an FCI, CF 
coordinates are the natural variables in terms of which 
the system can be understood in the simplest way. 



IV. THE CASE OF TRIVIAL BANDS 




Let us ask if we have achieved "too much" . Consider 
a topologically trivial band with Chern number zero. 
Nothing prevents us from representing its projected den- 
sity in terms of p G Mp(q + G) 



FIG. 9: The upper half shows the full landscape of the real 
part of the coefficient of d) ([p + q])d(p) in pFCB(q) in the 
approximation with 50 harmonics (from —25 to + 25) in the 
p y direction and 20 in the p x direction. The lower one shows 
the actual values. 



p(q) 



G 



c(G)p GMP (q+ G), 



(90) 



and carrying out the CF substitution. So, are there 
FCIs in topologically trivial bands? 

This appears to be a subtle issue. On the face of it, 
since the Bloch spinor can now be globally defined as a 
periodic function in the BZ, so can Pfcb, and its expan- 
sion in terms of p G mp which has a jump in the BZ seems 
ill fated. Although completeness assures us that with in- 
finite number of terms we can do it, in order to give the 
expansion the best chance, we must first transform the 
spinor to the gauge A y = 0, just like the functions enter- 
ing p GMP . This makes A x (p x ,p y + 2ir) ^ A x {p x ,p y ) and 
causes the familiar jump. 

The jump of this sort is inevitable if B(p) is non- 
constant, because in this gauge 

f p y 

■A x (p x ,p y ) = - B{p x ,p y )dp y (91) 
■'-f 

and this integral need not vanish at any fixed p x . 




FIG. 10: Comparison, in the case with C = 0, of imaginary 
parts of exact and approximate densities at p x = —5 + 1 
and p x = — If + 6 with the approximation truncated after 50 
harmonics (+25 to —25) in the p y direction 



The Fourier expansion, while not so successful as in 
the case C ^ 0, is still promising, as shown see Fig. [10] 
for two slices at p x = — | + 1 and p x = — | + 6. The 
key feature is that the smaller the jump (i.e. the smaller 
the magnitude of B(p)) the more Fourier components it 
takes to approximate it to a given accuracy. 

The bottom line is that, under certain conditions, even 
a band with C = could exhibit FQHE under partial 
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Hartree Fock Bands, nu=l/5 
Energy 



0.04 
0.02 

-0.02 
-0.04 




FIG. 11: The highest and lowest energies of the three sub- 
bands for the exotic fraction | as a function of the parameter 
Vpp. The bandwidths grows with the strength of the periodic 
potential but are still less than the sub-band gaps. The other 
parameters are V ee = .1,0) = 1. 



filling, an issue we are actively pursuing. 

So far we have ben able to rule out the following com- 
pletely trivial case: a band in which B(p) = with no 
band dispersion. The energy is all potential and propor- 
tional to V ee . Although an FCI state exists at v = |, its 
energy is always higher (by a factor of roughly 2) than 
the energy of the electronic Fermi liquid. Thus, in this 
case at least, FQHE does not appear to be favoured in 
the topologically trivial band. 

V. FRACTIONS WHERE v / a xy 

In the presence of a periodic potential, the equality 
v = G xy is not mandatory, since Galilean invariance is 
explicitly broken. Exotic states whose very existence 
depends on an interplay of interactions and the lattice 
potential are possible^. There is suggestive numerical 
evidence for such states in a pro blem of lattice hard-core 
bosons in an external field 35 36 . 

Since we have expressed the FCB Hamiltonian in CF 
language, it is natural for such states to be realized in 
FCBs under suitable conditions. 

Consider a case with v = \. Each electron sees 5 
flux quanta. Let us attach 2 (and not the usual 4) flux 
quanta to each electron to form CFs. Each CF sees 3 
flux quanta and v C¥ — | . In the continuum this partially 
filled CF-LLL is gapless, and thus not stable. However, 
as we will now show, such a state is generically gapped 
in the presence of a periodic potential. Since e*/e = 3/5, 
the CF sees | flux quanta in the electronic unit cell. So 
we must take 5 times the electronic cell (say in the x- 
direction) to form a CF unit cell that encloses integer 
flux (3) and ensures that the CF translation operators 
T x and T y commute. In general if there are p/q quanta 
of effective flux per unit cell, each CF-LL will split into 
p sub-bands^. Thus, in our example each CF-LL will 
spilt into three sub-bands. If the lowest of these three 
sub-bands, when filled, is separated from the others by a 
gap, we obtain an FCI. 

Our numerical work fully corroborates this picture. 
We work to linear order in the periodic potential Vpp, 



n= 1/5, V_ee=10,V_PP=.02 



s=0 




FIG. 12: From the top, the occupation numbers of the three 
CF-LLL- sub-bands, s=0,l,2, for the exotic fraction v = |. 
The action of T£ is to change p y by multiples of 

n? and change the sub-band index by 1 whenever p y changes 
by 2-7T. If the pictures of the three sub-bands are glued to 
together in the order 0, 1, 2, we see 5 oscillations altogether. 



and keep only the lowest CF-LL. Note that the treatment 
is not perturbative in the interaction strength. Fig. |11| 
shows the three sub-bands whose widths grow with the 
strength of the periodic potential while the gaps grow 
even faster. Figure [12] is very interesting. Once again 
T^T^T 3 ,^ commute with H but not T y . Each power 
of T x changes p v by Starting with p y = —tt of the 
sub-band s = 0, if we act with these powers of T x , we 
move to s — 1 when the momentum change is > 2-7T, and 
to s = 2 when the momentum change is > An. If the 
three pictures are glued end to end, we see 5 full oscil- 
lations. (Compare this to the v = | case where three 
periods occurred within the same CF-LL whereas here 
the 5 oscillations are spread over 3 sub-bands.) 
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What will be the Hall conductance of this state? We 
know from Kol and Read 3 that (in our convention where 
C = -1 for a filled LL) 



where cr CF , the dimensionless CF Hall conductance 
(equal to minus its Chern number in our convention) of 
the filled sub-band could be any integer. (When a LL 
splits into sub-bands we only know that the sum of the 
Chern numbers of the sub-bands equals that of the origi- 
nal LL. ) But no matter what this integer is, the possible 
values of <r xy = §,§,§■■ do not include |. The actual 
value happens to be | because a CF = 2. 

VI. SUMMARY AND OUTLOOK 

We have demonstrated here that the Hamiltonian the- 
ory of the FQHE 22 , which was very useful in describing 
the FQHE in the continuum, is just as effective in describ- 
ing fractionally filled Chern bands that exhibit FQH-like 
effects. This is surprising in view of the fact that since 
there is no external magnetic field in the FCB, ideas of 
flux attachment seem doomed from the outset. This is 
not an issue in our Hamiltonian theory which relies on 
an exact algebraic mapping that expresses the projected 
density of the Chern band p FCB (q) in terms of CFs in 
two steps. First, we express p FCB as a linear combina- 
tion of operators p GM p satisfying the magnetic transla- 
tion algebraic Second, we perform the CF substitution 
in Pgmp exactly as we did in the continuum theory 2 ^. 
These mappings of operator algebras can be carried out 
for arbitrary lattice models with no reference to Landau 
levels. 

We have illustrated the power of the approach by solv- 
ing two models. The first model is a v = | FCI and is adi- 
abatically connected to the continuum Laughlin state at 
the same filling. The second is a more exotic, v = | FCI 
relying on the lattice potential for its very existence^. 
Its dimensionless Hall conductance is | and not | as 
would be expected in Galilean invariant state. In both 
cases we have computed the band structure of CFs in the 
Hartree-Fock approximation. 



There are many interesting directions which we intend 
to pursue in future work. Collective excitations for frac- 
tionally filled Chern bands can be c omputed in a conserv- 
ing approximation in our approac h 1 22 ! 29 ! 30 !, as can finite 
temperature effects^!. 

A specially interesting case is at v = | in a FCB. One 
expects an electronic Fermi liquid at weak coupling and 
an HLR-type CF-Fermi liquid at strong coupling. This 
transition, which we propose to study by our operator- 
based method, has already been explored in the parton 
formulation recently^!. 

Let us turn to quasiparticle excitations. In addition to 
Laughlin-type quasiparticles with fractional charge and 
statistics, the lattice allows us to consider excitations not 
conceivable in the continuum, such as those associated 
with lattice vacancies or dislocations^. 

As stated in the introdu ction, fractionally filled 2D 
time-reversal invariant TIs^^^l can b e treated by la- 
belling the pair of Chern bands making up the nonintcr- 
acting TI (with equal and opposite Chern index) by a 
pseudospin index. There is no requirement of S z conser- 
vation, and the interactions can produce states which can 
spontaneously break time-reversal, and/or states of the 
Kol-Read type. In fact, in a numerical diagonalization on 
a small system Neupert et a¥^ find a state in a regime 
of parameters which has a filling of | but a degeneracy 
of 3 (rather than the degeneracy of 3 x 3 = 9 one would 
expect for "independent" v = | for each pseudospin) on 
the torus, suggesting that it could be a Kol-ReacP^ type 
state. 

While finishing this manuscript we noticed the work 
of Grushin et aF^I, who have examined the conditions for 
the stability of the FCI. 
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with the graphics. We are grateful to Shivaji Sondhi 
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